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Abstract Using the reference hypernetted chain (RHNC) integral equation theory and a rigorous stability analysis 
method, we investigate the phase behavior of a mixture of hard-sphere dipoles and neutral hard spheres based on 
the correlations of the homogeneous isotropic phase. Lowering the temperature down to the points where the RHNC 
equations fail to have a solution, several fluctuations strongly increase. At low densities our results indicate the onset of 
chain formation, which is similar with the pure DHS system. At high densities, the results indicate the appearance of 
isotropic-to-ferroelectric transitions (small neutral hard spheres concentrations) and demixing transitions (large neutral 


hard spheres concentrations). 
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Recently, several publications have discussed the sta- 
bility of fluids whose components carry permanent dipolar 
moments. It has been shown that the low temperature 
behavior of dipolar model fluids is much more complex 
than previously expected. For the most simple dipolar 
model fluids (i.e. dipolar hard sphere fluids), the sim- 
ulation results have shown that the system can sponta- 
neous polarize in dense liquid states!'!—*] and theoretical 
research has confirmed the appearance of a ferroelectric 
phase transition.!4—®] On the other hand, chain formation 
has been pronounced by recent theoretical research in the 
dilute gas states,[38-1°] the particles associate into chains 
and the conventional gas-liquid coexistence (which was as- 
sumed since the orientational interaction between two par- 
ticles is a van der Waals-like attraction!!!) is absent. 1” 
This behavior was also found by simulations.!3:9:18—15] 

The phase behavior of dipolar fluid mixtures is much 
more complex. As a simple model, DHS-HS mixture 
which has been considered by integral equation,{1®—18 
the results shown that the system demix within the 
isotropic phase at sufficiently large densities and coupling 
strengths. This phenomenon has been proved by Monte 
Carlo simulation.!!°! It seems that the effect of dipolar 
forces in mixture may be very different from that in pure 
fluids. Indeed, given that pure DHS fluids spontaneous 
polarize in dense liquid states and associated into chains 
with head-to-tail aligned dipoles in the dilute gas, an im- 
portant question is that whether DHS-HS mixtures can 
exhibit the similar phase behavior under certain condi- 
tions. 

To address these topics we employ reference hypernet- 
ted chain (RHNC) integral theory, which is used to cal- 
culate approximate two-particle correlation functions in 
the homogeneous isotropic fluid phase. Combining with a 
rigorous stability analysis which is on the base of diverg- 
ing fluctuations, we can characterize the phase transition 
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at the low temperatures. Previous application of this ap- 
proach have shown that the phase behavior is correctly 
predicted and the results are consist with the computer 
simulations.!”8:9] We note that the RHNC integral equa- 
tion theory is reliable to predict the low temperature phase 
behavior of DHS fluid mixtures. 

We consider a mixture of dipolar hard spheres (d) with 
dipolar moment u and neutral hard spheres (n), the two 
species have the same diameter ø. The pair potential be- 
tween two particles with coordinates (1) = (ri,w1) and 
(2) = (r2,w2) is given as 


o0 r2 <0, (1) 
uab(12) = ? 
Udip(r12, W1, w2); T12 >00, (1’) 

where rig = |ri2| = |r2 — rı| is the particle separation, 


w = (0, $) represents the orientation of a dipole in a spa- 
tially fixed coordinate system, and the subscripts a and 
b denote the components considered [a,b = n,d]. The 
dipolar potential is given by 


Udip (P12, #1, w2) = TA) . fiw) — 3[A(w1) - ĉ12] 


x [fwe) - ĉ12]}, (2) 
where ji(w) and fiz are unit vectors in the direction of w 
and ris. 

The full information about the system is contained 
in the total and direct correlation functions has(12) and 
Cab(12). We calculate these correlation functions in the 
homogeneous, high-temperature fluid phase via a numer- 
ical solution of the Ornstein-Zernike (O-Z) equationl?!! 


hap(12) = cas(12) +> f dBeae(13)pe(3)he(32), (3) 


combined with the RHNC closure approximation, 


Cap(12) = —1 — Nab (12), T12 S Cab, (4) 
Cab(12) = — ßbuab(12) = ln[1 + hab(12)] + hab(12) 
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+ Bav(12), T12 > Cab; (5) 


where Map(12) = hav(12) — can(12), 8 = (ks T)™! is the in- 
verse temperature (with ks being Boltzmann’s constant), 
and Bap(12) is so-called Bridge function, which are ap- 
proximated by the Verlet—Weis hard sphere correlation 
functions.” To solve equations (3)-(5) we expand the 
angle dependent functions in rotational invariants?’ 


fav(12) =y fle Pn www), (6) 


Iylol 


where f = h or c or u, O!2!(Wwow,) are the rota- 
tional invariants, a linear combination of the spherical 
harmonics. 23,24] 


DHe wwe) = XO Clll; mmm) 


mimm 


x D ow 


1)DRo(w2)Yim(wr). (7) 
The angels w1, w2, wp respectively describe orientations of 
the dipole at r1, r2, and the jointing vector r with respect 
to a space fixed coordinate system. The C(iylgl; miman) 
are the Clebsch—Gordan coefficients. Then finally after 
expansion, the RHNC equations combined with the O-Z 
relations become a system of integral equations for the r 
dependent coefficient, which is solved by iteration, here 
the expansion is truncated with 1,,l2 < 3. Details of the 
solution procedure can be found elsewhere. |!617:25—27] 

In the present work, our aim is to characterize the 
phase behavior of the system upon cooling from an 
isotropic high-temperature state. In order to do this, 
we apply a stability analysis developed by Chen and 
Forstmann!!”] which gives us a rigorous criterion and then 
we can infer the character of the phase transition. 

The central idea is to consider the fluctuation 6Q = 
Q-—O of the grand canonical free energy caused by small 
density fluctuation ôpa(1) around the homogeneous and 
isotropic equilibrium state, characterized by the equilib- 
rium density pa(1) = pa/4a. Up to the second order, 62 
is given by 


1 Q 
60 & oe È pr(2) 
= ks T pa 
dl i ie = can(12)] 
x Sal A (8) 


Stability (or metastability) implies that 6Q is positive for 
any small fluctuation with arbitrary direction in density 
space. 

Introducing Fourier transforms of the density fluctua- 
tions and correlation functions, respectively, and expand- 
ing these quantities along the k-frame (which is defined by 
choosing the polar axis along k), then Eq. (8) transforms 
to 


50 = EST oot IMC 


here the indices i = (a,1,) and j = (b, l2) is a combination 
of particle kind index and angular index, and the fluc- 
tuations dp;(x,k) = 5p!4-*(k)/,/pa- The elements of the 


k))ig D5 (x k) , (9) 
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matrices M[|x, k] are defined as 
M(x, k)i = 6:5 — (—1)*,/——~#-— 


x Gtx (k). (10) 


With the help of the expansion of the free energy fluc- 
tuation Eq. (9), we are able to study phase transitions 
like the gas-liquid condensation, demixing or orientational 
phase transitions in detail. The mean value of some fluc- 
tuations will diverge at the point of phase instability. But 
usually the condensation is accompanied by a change of 
concentration, or orientational phase transition and seg- 
regation appear simultaneously. Generally several fluctu- 
ations diverge at the phase transition. Only if dQ is a 
sum of pure squares of fluctuations, a vanishing prefactor 
of one fluctuation will lead to the divergence of just this 
fluctuation, which can be obtained by diagonalizing the 
matrix M in Eq. (9). The stable phase possess only posi- 
tive eigenvalues. The phase transition is indicated by one 
eigenvalue going to zero. The related eigenvector defines 
the combination of fluctuation, which will become singular 
at the transition point. 

For number density fluctuations we consider the 2 x 2 
submatrix Do of M(0,0) which is given as 


1- patna (0) —/PnPaend (0) 
Do = ( 0,0 ~00,0 ) : (11) 
-PERLO 1- pa (0) 

A vanishing of the smaller eigenvalue A, of Do indicates 
demixing or condensation of the system. This can be clar- 
ified by investigating the direction of the eigenvector re- 
lated to Ap. 

For polarization fluctuations, we consider the two 
eigenvalues Àr, of the 2 x 2 submatrix Fo of M(0,0) and 
Ar, of the matrix M (1,0). 


PaPb 
Ch 


J: 

Aro = 1 — >21. (0), (12) 
Pa 
Doug 

Ar, = 1+ =e)" (0) (13) 
Pd 


Ferroelectric transition of the system is indicated when 
eigenvalue Àr, goes to zero. To see this we consider the 
physical precursor of ferroelectric order — that is, the mix- 
ture’s dielectric constant e. The relation of this quantity 
and the two eigenvalues isl8—18] 


Ne 
es 


Numerical investigation shows that Ar, remains finite at 
all states considered. Therefore, a vanishing of Ar, can be 
uniquely related to a divergence of the dielectric constant. 


(15) 
which indicates an isotropic-to-ferroelectric transition. 
The state of the mixtures of dipolar hard spheres and 
neutral hard spheres is characterized by three physical 
parameters: the reduced dipole moment u* defined by 
pr? = p*/o3kgT = 1/T* (T* is the reduced temper- 
ature), the reduced density of the whole mixture p* = 
(Pn + pa): o° = p- o? and the concentration of neutral 


E€ = 


(14) 


e> oœ & Ap, — 0, 


No. 3 


particles c = Cn = pn/p. Lowing the temperature at con- 
stant total density and concentration we find a tempera- 
ture T3(p*, cn) below which the RHNC equation have no 
solution, where the homogeneous isotropic mixtures be- 
come unstable due to a phase transition. In the following 
we discuss results obtained by varying the concentration 
of neutral particles and temperature at two representative 
values of p*, one is dilute system with very low densities 
p* = 0.04, and the other is hight densities p* = 0.8. 

We start by considering the mixtures at p* = 0.04. 
The behavior of pure dipolar hard sphere fluids with 
p* = 0.04 is investigated by the RHNC integral equation 
approach!! and the Monte Carlo simulations.!3-9! The re- 
sults from these two methods all indicate that at low den- 
sities and certain low temperatures dipolar hard spheres 
form pairs and chainlike clusters instead of an ordinary 
gas-liquid transition. 

Turning now to mixtures of dipolar hard spheres and 
neutral hard spheres, we determine the stability limits of 
the mixed isotropic phase — i.e., the temperatures T* (cn) 
where the RHNC equations have no solution lowering T* 
towards Tž. Here we consider two representative concen- 
trations cn = 0.1 and cn = 0.8. These two systems con- 
sidered indeed exhibit pair formation which can be seen 
most directly via the pair correlation functions h2%?(r) 
(see Figs. 1 and 2). These functions are the projections of 
the total correlation functions hgg(12) onto the spherical 
invariant ®?°? defined by 


1 

02 2 TBG) fra)? 1. 
From the definition of ©°, large values of h30? (r) thus 
indicate that the dipolar hard spheres in the neighbor- 
hood of a certain dipolar hard sphere arrange preferably 
at the “poles” of the sphere. Both in Figs. 1 and 2 we 
see that the first maximum in h20? (r), located at r = ø, 
grows drastically for T* — Tž, it reflects the growing ten- 
dency of two dipolar hard spheres in contact to have a 
“nose-to-tail-alignment”. Such an arrangement is not sur- 
prising because it represents the configuration with the 
lowest energy of the system. 

We now compare our results with those for the pure 
dipolar hard sphere system (cn = 0) from RHNC approach 
by S. Klapp and F. Forstmann.!l In the case p* = 0.04, 
Cn = 0, the value of T is 0.329, which is higher than the 
case Cn, = 0.1 and cn = 0.8, we have Tž (cn = 0) > TF 
(Cn = 0.1) > Ti (cn = 0.8). We understand this as a con- 
sequence of increasing frustration effects: the more dipo- 
lar hard spheres are replaced by the neutral hard spheres, 
the more the two particle interaction which leads to pair- 
ing and cluster formation, becomes disturbed. Therefore 
the temperatures which are necessary to reach the “chain- 
state” become lower. 

To discuss the phase behavior of mixtures at high den- 
sities we consider the value p* = 0.8 as a representative 
example. We first determined the instability line of the 
systems at p* = 0.8 for each concentration c,. Results 
are shown in Fig. 3. Upon lowering T* towards T, the 
RHNC equation have no solution. For the case cn = 0, the 
system is again pure DHS fluids. Turning back to Fig. 3, 


(16) 
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our results show that the ferroelectric instability appears 
on the left-hand side of the diagram (small cn), just as in 
the pure DHS system. But on the right-hand side of the 
diagram (large cn), the isotropic demixing appears. This 
is evidenced by the inverse temperature dependence of the 
eigenvalues Àr, [see Eq. (21)] and Ap [see Eq. (17)], which 
are plotted in Fig. 4 for the representative case of cn = 0.2 
and cn = 0.8. 
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Fig. 1 The coefficient h28? (r) at p* = 0.04, cn = 0.1, 
and T* — Tg. The inset shows the development of the 
maximum of the first peak vs inverse temperature. 
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Fig. 2 The coefficient h2% (r) at p* = 0.04, cn = 0.8 
and T* — Tg. The inset shows the development of the 
maximum of the first peak vs inverse temperature. 


The results in Fig. 4(a) clearly show that Ap, goes 
to zero significantly faster than Ap at the case cn = 0.2, 
the polarization fluctuations become dominating for T* — 
Tg. The system thus transform directly into ferroelec- 
tric phases as do in the pure DHS system. On the other 
hand, at the large concentration cn = 0.8 (Fig. 4(b)) Ap 
goes to zero firstly, and using the direction of the related 
eigenvector , indicating that the diverging fluctuations re- 
lated to Ap are essentially concentration (rather than total 
density) fluctuations. All our results are consistant with 
the RHNC study of monodisperse DHS mixtures!!§) with 
small interaction parameters. 

Finally, it is worth to note from Fig. 3 that on the left- 
hand side of the diagram, the instability temperatures Tè 
are slowly increased upon replacing dipolar hard spheres 
more and more by neutral hard spheres. It indicates that 
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the appearance of the neutral hard spheres is helpful for 
a ferroelectric transition. 
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Fig. 3 The instability line of dipolar and neutral hard 
spheres mixtures with the total density p* = 0.8 in the 
concentration-temperature plane. 


In conclusion, we have seen that for DHS and neutral 
hard sphere mixtures, the RHNC theory is capable to pre- 
dict the phase transition both at small and high densities. 
At very low total density p* = 0.04 which is in the range 
where the pure DHS system associate into short clusters 
at sufficiently low temperatures. The same type of phase 
behavior is found to occur in our DHS and neutral hard 
sphere mixtures for all concentration c,, and the insta- 
bility temperature T% is gradually decreased when more 
and more dipolar hard spheres are replaced by the neu- 
tral hard spheres. At the large total density p* = 0.8, our 
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